/*
  Copyright (C) 2011, Kyungjoo Kim
  
  This file is part of UHM (Unassembled HyperMatrix) sparse direct solver.
  
  All rights reserved.

  This library is free software; you can redistribute it and/or
  modify it under the terms of the GNU Lesser General Public
  License as published by the Free Software Foundation; either
  version 2.1 of the License, or (at your option) any later version.
  
  This library is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  Lesser General Public License for more details.
  
  You should have received a copy of the GNU Lesser General Public
  License along with this library; if not, write to the Free Software
  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
  Also add information on how to contact you by electronic and paper mail.
  
  
  Kyungjoo Kim
  iamkyungjoo@gmail.com
*/
#ifndef UHM_WRAPPER_FORT_LINAL_DOT_PETSC_HXX_
#define UHM_WRAPPER_FORT_LINAL_DOT_PETSC_HXX_


extern "C" {

  // ** Solvers
  extern F_void UHM_C2F(uhm_petsc_lu_nopiv_initialize)(F_int *type, F_int *rhs, F_int *blk,
                                                       F_ptr *solver);
  extern F_void UHM_C2F(uhm_petsc_lu_piv_initialize)(F_int *type, F_int *rhs, F_int *blk,
                                                     F_ptr *solver);
  extern F_void UHM_C2F(uhm_petsc_lu_incpiv_initialize)(F_int *type, F_int *rhs, F_int *blk,
                                                        F_ptr *solver);
  
  extern F_void UHM_C2F(uhm_petsc_ldl_nopiv_initialize)(F_int *type, F_int *rhs, F_int *blk,
                                                        F_ptr *solver);
  extern F_void UHM_C2F(uhm_petsc_ldl_piv_initialize)(F_int *type, F_int *rhs, F_int *blk,
                                                      F_ptr *solver);
  extern F_void UHM_C2F(uhm_petsc_ldl_incpiv_initialize)(F_int *type, F_int *rhs, F_int *blk,
                                                         F_ptr *solver);

  extern F_void UHM_C2F(uhm_petsc_chol_nopiv_initialize)(F_int *type, F_int *rhs, F_int *blk,
                                                         F_ptr *solver);
  extern F_void UHM_C2F(uhm_petsc_chol_piv_initialize)(F_int *type, F_int *rhs, F_int *blk,
                                                       F_ptr *solver);
  extern F_void UHM_C2F(uhm_petsc_chol_incpiv_initialize)(F_int *type, F_int *rhs, F_int *blk,
                                                          F_ptr *solver);

  // ** PETSC Setting
  extern F_void UHM_C2F(uhm_petsc_assemble)(F_ptr *solver);
  extern F_void UHM_C2F(uhm_petsc_set_alloc_scale)(F_ptr *solver, F_int *alloc);
  extern F_void UHM_C2F(uhm_petsc_set_tolerence)(F_ptr *solver, F_double *rtol, F_double *atol);
  extern F_void UHM_C2F(uhm_petsc_set_pctype)(F_ptr *solver, F_char *pctype);

}


#endif
